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Abstract. 

Nearby clusters and groups of galaxies are potentially bright sources of high-energy 
gamma-ray emission resulting from the pair-annihilation of dark matter particles. 
However, no significant gamma-ray emission has been detected so far from clusters 
in the first 11 months of observations with the Fermi Large Area Telescope. We 
interpret this non-detection in terms of constraints on dark matter particle properties. 
In particular for leptonic annihilation final states and particle masses greater than 
~ 200 GeV, gamma-ray emission from inverse Compton scattering of CMB photons 
is expected to dominate the dark matter annihilation signal from clusters, and 
our gamma-ray limits exclude large regions of the parameter space that would 
give a good fit to the recent anomalous Pamela and Fermi-LAT electron-positron 
measurements. We also present constraints on the annihilation of more standard dark 
matter candidates, such as the lightest neutralino of supersymmetric models. The 
constraints are particularly strong when including the fact that clusters are known 
to contain substructure at least on galaxy scales, increasing the expected gamma-ray 
flux by a factor of ~ 5 over a smooth-halo assumption. We also explore the effect of 
uncertainties in cluster dark matter density profiles, finding a systematic uncertainty 
in the constraints of roughly a factor of two, but similar overall conclusions. In this 
work, we focus on deriving limits on dark matter models; a more general consideration 
of the Fermi-LAT data on clusters and clusters as gamma-ray sources is forthcoming. 
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1. Introduction 

Clusters of galaxies are the most massive objects in the Universe that have had time 
to virialize by the present epoch, making nearby clusters attractive targets for searches 
for a signature from dark matter annihilation. Clusters are more distant, but more 
massive than dwarf spheroidal galaxies, and like dwarf spheroidals, they are very dark 
matter dominated and typically lie at high galactic latitudes where the contamination 
from Galactic gamma-ray background emission is low. A common feature of the stable 
particle products resulting from the annihilation of weakly interacting massive particles 
(WIMPs), prime candidates for dark matter, is a significant gamma-ray emission, though 
in addition dark matter annihilation is expected to yield a broad multiwavelength 
spectrum of emission [I]. For a variety of particle models, the gamma-ray signal from 
the annihilation of dark matter in clusters was predicted to be potentially observable 
by the Fermi Gamma- Ray Space Telescope [21 El- 

In proto-typical WIMP models, such as for neutralinos in the minimal 
supersymmetric extension to the Standard Model, gamma-ray emission results primarily 
from the decay of neutral pions produced as part of the annihilation products. For 
WIMPs annihilating primarily into leptons, such as the dark matter models proposed to 
explain the Pamela positron excess [Jj, gamma-ray emission results instead dominantly 
from internal bremsstrahlung in the final state. In addition, these models lead to a large 
population of energetic e + and e~ among the annihilation products (hence the popularity 
as an explanation for the recent e + e~ measurements) which in turn lead to significant 
gamma-ray emission through the inverse Compton scattering (IC) of cosmic microwave 
background photons (CMB) or other photon fields in the inter-galactic medium. The 
IC emission is particularly important in clusters which effectively confine cosmic ray 
electrons much longer than the time it takes them to lose energy through this mechanism. 
We will show that Fermi-LAT observations of clusters place some of the best constraints 
on the possible interpretation of the local e + e~ signal observed by Pamela, Fermi-LAT, 
and H.E.S.S. [4J, [6] as a signal from the annihilation of dark matter [7J [8j [9]. 

Dark matter annihilation is not the only potential source of gamma-ray emission 
from clusters of galaxies. Clusters are also expected to host a population of cosmic rays 
accelerated in merger and accretion shocks and by AGN throughout the cluster's history. 
Collisions of relativistic protons with protons in the thermal intracluster medium would 
lead to gamma-ray emission through the decays of neutral pions produced in the 
collision. In addition, diffuse radio emission observed in many clusters reveals the 
presence of relativistic electrons, and similar to the dark matter case, IC scattering 
of background radiation by these relativistic electrons could also lead to detectable 
gamma-ray emission. The shape of the gamma-ray spectrum is generally expected to 
be different than that from dark matter annihilation, although sufficient statistics are 
required to distinguish between the two [2]. 

Significant gamma-ray emission has not been detected from local clusters by Fermi- 
LAT in the first 11 months of observations [TO], aside from the detection of the central 
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radio galaxies in a couple of clusters [TTJ [T2]. The Fermi upper limits on gamma- 
rays from clusters place interesting limits on dark matter models. In this paper, we 
investigate in detail the implications of the Fermi-LAT data in terms of limits on some 
of the fundamental properties of dark matter, such as the particle mass, pair-annihilation 
final state and pair- annihilation rate. We use 11 months of Fermi-LAT survey mode data 
and self-consistently derive the flux upper limits for dark matter models based on the 
expected gamma-ray spectrum as a function of particle mass, annihilation final state, 
and emission mechanism (Sec. [2]). For specific models, we present the Fermi exclusion 
regions for the best-candidate clusters in the annihilation cross-section - mass plane 
(Sec. [3]). In this work, we focus specifically on clusters as potential sources of a dark 
matter signal. The Fermi-LAT results for a much larger sample of clusters as well as a 
more general consideration of clusters as gamma-ray sources, including constraints on 
their cosmic ray content, will be presented in a forthcoming paper. 

2. Data Analysis 

2.1. Cluster Sample 

We selected our cluster sample from among the clusters predicted to have the brightest 
gamma-ray emission from dark matter annihilation based on the modeling described in 
[2]. The predicted gamma-ray flux from dark matter annihilation $ 7 is proportional to 
the dark matter density squared integrated along the line-of-sight and averaged on a 
given angular scale: 



In turn, for a distant object whose spatial extent is smaller than the angular scale under 
consideration, one has a simple scaling with the cluster distance (even though in the 
present study we resort to the full treatment in the equation above): 



where Vol indicates the volume of the cluster, and D is the cluster distance. In 
determining the brightest predicted clusters, we used the HIFLUGCS [13] catalog of 
the brightest X-ray clusters in the sky and X-ray determinations of the cluster masses. 
Starting from the best candidate clusters, we removed from the sample the Perseus 
cluster since its central galaxy, NGC 1275, hosts a bright gamma-ray AGN [11], the 
Ophiuchus cluster, which lies near the Galactic Center, and the Norma cluster, which 
lies near the Galactic plane. Our final sample of six clusters is listed in Table [U 

2.2. Fermi-LAT Observations and Reduction 

The Fermi-LAT data selection is the same as described in [15]. In brief, we utilize 
the first 11 months of scanning mode observations collected between August 4, 2008 
and July 4, 2009. Fermi operates in sky survey mode, covering the entire sky every 




(1) 




(2) 
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Table 1. Clusters considered in the Fermi-LAT analysis. Cluster positions 
are taken from [13], and redshifts from the NASA/IPAC Extragalactic Database 
(http: /7nedwww.ipac.caltech.edu/| . The Fermi-LAT data is analyzed within a 10° 
radius region surrounding each cluster position. Column 5 lists J, the integral along 
the line of sight of the cluster dark matter density squared, for an assumed NFW profile 
modeled as described in sec. l3.1l and including the smooth halo only (i.e. neglecting the 
potential effect of substructure). J is calculated within an radius of one degree (i.e. 
within a solid angle region Ail corresponding to one degree). Centaurus and NGC 



4636 are cool core clusters in X-ray [14] . 


Cluster 


RA 


Dec. 


z 


J (10 17 GeV 2 cm- 5 ) 


AWM 7 


43.6229 


41.5781 


0.0172 




Fornax 


54.6686 


-35.3103 


0.0046 


6-8±fcS 


M49 


187.4437 


7.9956 


0.0033 


4-4+0;? 


NGC 4636 


190.7084 


2.6880 


0.0031 


4 1+0-3 


Centaurus (A3526) 


192.1995 


-41.3087 


0.0114 


9 7+0.I 
z - ' -0.1 


Coma 


194.9468 


27.9388 


0.0231 


i Y+o.i 
L - 1 -0.1 



3 hours and providing nearly uniform coverage. We employ the standard data cuts 
and background rejection defined in [16] keeping only "Diffuse" class events which have 
a high probability of being gamma-rays. In addition, to avoid contamination from the 
Earth's "albedo" gamma-ray emission.)., we exclude time intervals where the zenith angle 
was greater than 105° (the Earth's limb lies at a zenith angle of 113° for Fermi) and 
times when the rocking angle (angle between the viewing direction and the zenith) was 
greater than 43°. The data were further limited to include only events with energies 
between 100 MeV and 100 GeV, where the LAT instrument is best calibrated and where 
models of the diffuse gamma-ray backgrounds are available^ and only data within a 10 
degree radius of each cluster. 

The data analysis was carried out using the Fermi Science Tools package version 
9.15.2 and using the Pass 6 version 3 LAT instrumental response functions (IRFs). 
The LAT data were fit using an unbinned likelihood method p21 [18] as implemented 
in the pyLikelihood task in the Science Tools. This task fits the LAT data to a 
combined spatial and spectral model of the source folded with the LAT response and 
point spread function (PSF) using a maximum likelihood statistic. To account for the 
gamma-ray backgrounds, the fit includes models for both the Galactic and isotropic 
diffuse gamma-ray emission^ gllJem_v02.fit and isotropic_iem_v02.txt, respectively) as 
well as surrounding point sources within a 15 degree radius. The diffuse emission 
models are those currently advocated and released by the Fermi-LAT collaboration, 
and the normalizations of these components are allowed to vary during the fitting. The 
background point sources are drawn from the First Fermi-LAT source catalog, and the 

| gamma-ray emission from cosmic ray interactions with the Earth's atmosphere 

§ http://fermi.gsfc.nasa.gov/ssc/data/analysis/LAT_caveats.html 

|| http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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parameters of these sources are left fixed at their best-fit catalog values [19^ 1 

No significant detections of gamma-ray emission from clusters of galaxies have been 
found in the first 11 months of Fermi-LAT observations [10] as will also be discussed in 
a forthcoming paper. In general, the background model including the Galactic diffuse, 
isotropic diffuse, and nearby point sources gives a good fit to the data. The most 
noticeable residual in the fits is the tendency for the model to slightly overpredict the 
data by around 10% at the lowest energies considered (100-150 MeV), likely due to the 
rapidly changing effective area at these energies. However, for the models discussed in 
this paper the flux upper limits we derive are essentially unaffected. The improvement 
in the LAT sensitivity and PSF at higher energies and the softness of the diffuse 
backgrounds mean that the dark matter models we consider are primarily constrained 
by the LAT data at higher energies. 

In order to place upper limits on the possible gamma-ray flux from dark matter 
annihilation, we model each cluster as a Fermi point source with an assumed dark matter 
annihilation spectrum. In particular, we fit the data for a grid of assumed dark matter 
particle masses and a couple of representative annihilation final states. Upper limits on 
the flux are calculated using a profile likelihood technique: values of the flux are scanned 
fitting with respect to the other free parameters and the change in the log-likelihood is 
calculated. The 95% confidence upper limit on the flux is defined as the flux value for 
a change in the log-likelihood of 1.36. 

2. 3. Flux upper limits for selected dark matter models 

We analyze the Fermi-LAT data to derive upper limits on the gamma-ray flux from 
the clusters of galaxies listed in Table 1. The Fermi flux limits depend on the shape of 
the assumed gamma-ray spectrum, and we consider dark matter annihilation in specific 
representative models, self-consistently incorporating the expected spectral shape for a 
given particle mass and annihilation final state. In Sec. [3], we will describe the translation 
of the gamma-ray flux limits to limits on the dark matter annihilation cross-section 
including modeling of the cluster dark matter density profiles. 

Fig. [T]shows the 95% confidence upper limits on the flux in the 100 MeV to 100 GeV 
band for the gamma ray spectrum produced by pair-annihilation into a bb final state 
versus particle mass. A bb final state is chosen as a prototypical case for supersymmetric 
dark matter. Gamma-ray emission results primarily from the decay of neutral pions 
among the annihilation products, and the expected gamma-ray spectrum is very similar 
for other quark-antiquark and gauge boson final states. Here the spectra are modeled 
using the DMFIT package [20J which has been incorporated in the likelihood packages 
in the Science Tools. As shown in Fig. [H the flux upper limits decrease with increasing 
particle mass. For larger particle masses and a bb final state, the expected dark matter 
annihilation spectrum in the Fermi-LAT energy range becomes significantly harder. 
The diffuse backgrounds are relatively soft and at higher energies the LAT PSF and 

http://fermi.gsfc.nasa.gov/ssc/data/access/lat/lyrj3atal0g/ 
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Figure 1. Flux upper limits as a function of particle mass for an assumed bb final 
state for the clusters listed in Table [T] The upper limits reflect the 95% confidence 
level limits on the photon flux in the 100 MeV to 100 GeV energy range. 



sensitivity improve, leading to stronger limits on models which predict a large number 
of high energy photons. The flux limits vary a bit from cluster to cluster as might be 
expected given variations in the diffuse backgrounds across the sky as well as variations 
in the proximity of bright gamma-ray sources. 

We also consider dark matter annihilating to a final state as an example of 

models giving a good fit to the Pamela anomalous positron fraction data [3] and to 
the Fermi-LAT e + e~ data [5]. (We do not explicitly consider here the case of a t + t~ 
final state, for which we would have obtained intermediate limits between the 
and bb cases.) For a leptonic final state such as gamma-ray emission results from 

both internal bremsstrahlung emission from the charged final state, referred to as final 
state radiation (FSR), as well as from Inverse Compton scattering (IC) of background 
radiation fields off of the energetic e + and e~ produced among the annihilation products. 
Here we include both sources of emission, but for the IC component, we conservatively 
only include IC scattering of the CMB. Scattering of other background radiation fields 
(such as starlight) could also contribute but are not expected to dominate for these dark 
matter models. In clusters, energetic e + and e~ are expected to lose energy, primarily 
through IC scattering, much more quickly (energy loss time scales of order 10 14 sees or 
less for E e ± > 100 GeV) than the time it takes for them to diffuse out of the system 
with any reasonable cosmic ray diffusion setup. The effect of cosmic rays escaping the 
diffusion region can thus effectively be neglected for clusters, as shown e.g. in Ref. pp. 

A suppression of the IC signal might occur if synchrotron losses, averaged over 
the entire volume of the cluster, are comparable or more significant than IC losses. 
Defining an effective average cluster magnetic field B, this condition amounts to 
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Figure 2. Left: The expected gamma-ray spectrum from dark matter annihilation 
into a /J + /i~ final state including both FSR and emission from IC scattering of CMB 
photons for different assumed particle masses. The spectral normalization is arbitrary, 
but the same effective J is used in all cases. Right: Flux upper limits as a function 
of particle mass for an assumed final state, including the contributions of both 

FSR and IC gamma-ray emission, for the clusters listed in Table [TJ The upper limits 
reflect the 95% confidence level limits on the photon flux in the 100 MeV to 100 GeV 
energy range. 



B > -Bcmb — 3.2(1 + z) 2 /zG, the latter quantity being the amplitude of a magnetic 
field having the same energy density as the CMB. While such large magnetic fields are 
possible, they are unlikely to uniformly populate the cluster volume we consider. If, 
however, this were the case, the suppression in the IC signal we would get (and thus 
the relaxation of the constraints we present in Sec. [3] below) would be on the order of 

In the Fermi-LAT energy range, the IC emission is expected to dominate for dark 
matter particles annihilating into n + fi~ for particle masses greater than a couple hundred 
GeV. As DMFIT does not include the contribution of IC, we fit the Fermi-LAT data 
using custom spectra for different assumed particle masses including both the IC and 
FSR gamma-ray radiation. Fig. |2] (left) shows examples of the expected spectra for 
different particle masses. Fig. [2] (right) shows the corresponding upper limits on the 
gamma-ray flux as a function of particle mass. For the lowest masses (less than ~ 200 
GeV) and the highest masses (greater than ~ 1 TeV) the expected spectrum is fairly 
hard; at low masses FSR dominates the emission with a typical dN/dE ~ E^ 1 spectrum 
where at high masses the IC spectrum itself is fairly hard. The largest upper limits come 
for particle masses of a few hundred GeV where the IC emission begins to dominate but 
the spectrum is fairly soft. 

In calculating the upper limits, we have assumed that the clusters are point sources 
for Fermi. The gamma-ray signal from dark matter annihilation is proportional to the 
dark matter density squared making it centrally concentrated. In the next section, we 
will describe in detail the modeling of the cluster dark matter density distributions for 
which we assume an NFW profile. For the clusters in our sample, which have NFW 
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scale radii between 0.25 and 0.45 degrees, the point source approximation is reasonable. 
However, particularly for very hard spectra which predict a significant number of high 
energy photons where the Fermi-LAT PSF is smaller, some clusters could be slightly 
extended for Fermi. We test the effect of modeling clusters with an NFW spatial model 
compared to a point source model for a hard power law spectrum with an index of 1. 
We find that our upper limits are only expected to increase by 10 — 40% when taking 
in to account the spatial distribution, which is less than the effect of the uncertainty 
in the cluster density profiles discussed in Sec. 13.31 Another possible source of error 
in the fluxes comes from the systematic uncertainty in the LAT effective area. This 
uncertainty is energy dependent and has been estimated to be 10% at 100 MeV, 5% at 
560 MeV, and 20% at 10 GeV [22]. 



3. Results 



3.1. Constraints on Dark Matter Annihilation Models 

In this section, we explore the constraints on models of dark matter annihilation implied 
by the non-detection of nearby clusters of galaxies with Fermi-LAT. The expected 
gamma-ray flux from dark matter annihilation in a given object depends on both a 
particle physics factor, for a given annihilation final state or combination of final states 
the normalization is proportional to the annihilation cross-section divided by the particle 
mass squared, and an astrophysical factor dependent on the integral of the dark matter 
density distribution squared along the line of sight, J given in Eq. 1. 

Following the results of cosmological collisionless N-body simulations, we assume 
a Navarro- Frenk- White profile [23] for the dark matter density distribution p(r) of the 
clusters in our sample, where r is the radial distance from the center of the cluster 
(which we assume to be spherically symmetric): 

PNFw( r ) = — f (3) 

(r/r s ) (1 + r/r s ) 

We calculate the values of the scaling density p s and scaling radius r s from the 
cluster virial mass, as determined from X-ray observations in Ref. [13], and from the 
observationally determined concentration-mass relation of Ref . [24J: 

where the virial concentration c v i r = Rvn/^s for an NFW profile. The J values, within 
an angular region of one degree radius, determined by this method are listed in Table [TJ 
The errors on J listed in Table [1] reflect the uncertainties in the X-ray virial mass 
determinations [13], but below we will also discuss alternate determinations of the 
density profiles for a couple of the best candidate clusters. 

Fig. El shows the upper limits on the annihilation cross-section for a bb final state 
(left panel) and a fi + fi~ final state (right panel) as a function of WIMP mass for all 
clusters given the upper limits on their gamma-ray fluxes (Fig. [T] and Fig. |2]) and the 
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Figure 3. Upper limits on the annihilation cross-section for a bb final state (left panel) 
and a final state (right panel) as a function of WIMP mass for all clusters in 

the sample. In the left panel, the constraints are compared to the predicted cross 
section for generic SUSY WIMP models (yellow) in the general MSSM [25] as well as 
to the subset of models that have thermal relic densities compatible with the observed 
universal dark matter density (orange) |26j . In the right panel, the cluster gamma-ray 
constraints are compared to values of the mass and annihilation cross section for dark 
matter models annihilating to a \i + \x~ final state which provide good fits to the Pamela 
(blue) or Fermi-LAT (orange) e + e~ data [7]. 



modeling of their dark matter density profiles. The tightest constraints are given by 
nearby groups of galaxies, Fornax, M49, and NGC 4636, which have the largest J 
values. For the case of a bb final state with no contribution from substructure, the 
cluster constraints exclude only models with relatively high annihilations cross- sect ions 
(which would produce low thermal relic densities) and are generally weaker than those 
found for nearby dwarf spheroidal galaxies (dSph) [T5] . 

For a final state, models fitting the Pamela positron excess with masses above 

2 — 3 TeV are ruled out by the non-detection of the Fornax group. Similar constraints 
were derived using the Ursa Minor dSph galaxy [15] for a diffusion coefficient of D = 10 28 
cm 2 /s. For a larger effective diffusion coefficient, cosmic rays would leave the dwarf 
system at much larger rates, and the IC emission would be significantly reduced, giving 
weaker constraints for dSphs compared to clusters. As argued in Sec. 2.3, the constraints 
shown here for clusters are therefore more robust in that the cosmic ray escape rate for 
high-energy e + and e~ is virtually negligible for cluster size systems. 

Our results improve significantly the constraints on particle dark matter obtained 
from observations of nearby clusters of galaxies by Atmospheric Cherenkov Telescopes, 
namely observations of the Coma cluster by H.E.S.S. [27] and of the Perseus cluster 
by MAGIC [28]. Notice that we do not consider here the search for monochromatic 
gamma-ray lines, since the predicted constraints from observations of galaxy clusters 
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are much weaker than those obtained in the dedicated search of Ref . [2H] • 

3. 2. Effects of Substructure 

The above constraints all assumed only the smooth cluster halo contribution to the 
dark matter density. However, N-body simulations in a CDM cosmology show a large 
spectrum of dark matter substructures in cluster and galaxy halos [301 El] • F° r typical 
thermally produced WIMP models, these substructures are expected to exist down to 
cutoff scales of between 10~ 6 and 10 2 Earth masses [32j [331 El]- These dense dark 
matter substructures can add significantly to the gamma-ray flux from annihilation. At 
the very least, massive groups and clusters are known to host hundreds to thousands of 
galaxies ranging in size from giant ellipticals down to dwarf galaxies. 

We estimate the impact of substructures on J (the integral in Eq. (T2])) and the 
corresponding boost in the expected gamma-ray flux from clusters in our sample using 
the following procedure. First, we assume the concentration mass relation for small scale 
halos found in the numerical simulations of Ref. [35] , which contains a radial dependence 
accounting for the tidal stripping of substructure close to the center of the host halo. 
We assume the radial scale of said dependence to be twice the virial radius of the host 
halo. We also verified that neglecting the radial dependence of the concentration mass 
relation for subhalos only affects the resulting substructure enhancement by less than 
20%. We assume a spectral index of —1.9 for the subhalo mass function, and we consider 
two realizations for the substructure distribution: 

• A conservative setup where the fraction f s of cluster mass in substructure equals 
10%, and where we assume no substructure (i.e. we effectively cut-off the integration 
on the subhalo mass function) below the mass scale of a dwarf galaxy dark matter 
halo, M cut = 10 7 Mq. We consider this setup to be conservative as we know that 
clusters at least contain galaxies. 

• An optimistic setup, where f s = 0.2 and M cut = 10~ 6 M [36] (even smaller cutoff 
scales are however possible in well-motivated WIMP scenarios [32] - see also [31] 
for a review). 

These substructure setups in terms of the mass function of subhalos, the mass fraction in 
subhalos, and the subhalo mass-concentration relation were chosen based on the results 
of recent N-body simulations [311 EDI EZ]- We verified that the substructure setups 
we consider here are consistent with the alternative framework presented in Ref. [I]. 
Specifically, the optimistic setup corresponds to a substructure contrast factor A 2 (in 
the notation of Ref. [1]) of A 2 ~ 3 x 10 5 , while the conservative setup to one with 
A 2 ~ 1.3 x 10 5 , for the given specified substructure mass fractions. 

Fig. H] shows the effect of adding substructure to the dark matter constraints shown 
in Fig. [3] for two example clusters, Fornax (which has the highest J for a smooth halo) 
and Coma (the most massive cluster in the sample). As in the previous section, we show 
the constraints for a pure bb final state (left panel) and a final state (right panel). 

The dashed lines show the results for our conservative substructure setup which includes 
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Figure 4. Upper limits on the annihilation cross-section for a bb final state (left panel) 
and a [i + n~ final state (right panel) for the Coma and Fornax clusters including the 
effect of substructure on the expected gamma-ray signal. The constraints are shown 
for no substructure (solid lines), our conservative substructure setup which includes 
only substructure of dwarf galaxy mass or larger (dashed lines), and our optimistic 
setup which includes substructure down to M cut = 10~ 6 M Q (dot-dashed lines). The 
dark matter models are the same as in Fig. [3] 



only the expected contribution from galactic scale substructure and gives a boost to the 
expected gamma-ray flux of ~ 4.6 for Fornax and ~ 2.1 for Coma. In this setup, the 
Fornax constraints exclude models fitting the Pamela e + e~ data with masses above 1 
TeV for a fi + fi~ final state and begin to probe thermally produced WIMP models with 
a relic density consistent with the observed universal matter density for a bb final state. 

In the more optimistic substructure setup (dot-dashed lines in Fig. HI), where we 
include substructure down to roughly the expected substructure cutoff scale for WIMP 
models, the predicted boosts to the cluster gamma-ray flux are even higher, ~ 10 for 
Fornax and ~ 9 for Coma. In this case, the left panel indicates that winos lighter than 
200 GeV are ruled out, as is a dark matter interpretation of the Pamela positron fraction 
with a fi + n~ final state and a sub- TeV mass. In general, models with a final state 

with cross sections greater than ~ 10~ 23 cm 3 s _1 and particle masses greater than ~ 1 
TeV are excluded which includes most of the parameter space fitting the Fermi-LAT 
e + e~ data. 

Similar to our conclusions here, Ref. [3] predicts that for a particle mass of 1.6 
TeV with a cross-section of 3 x 10 -23 cm 3 s^ 1 annihilating to a pure final state, 

Fermi-LAT should detect local clusters if the cut-off scale for substructures is 10 4 M Q or 
less. 
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3.3. Uncertainty in the Density Profiles 

In this section, we consider alternative determinations of cluster dark matter density 
profiles. Our primary results are based on cluster dark matter density profiles from X-ray 
observations, and Table [T]lists the uncertainties in J stemming from the X-ray virial mass 
determinations [13J only, typically about 10% in J. However, many of the clusters in 
the sample have masses and density profiles determined through other observables, and 
these perhaps better reflect the total uncertainties. For example, the density distribution 
of the Coma cluster has been determined through both weak lensing [38, 39] and galaxy 
dynamics [ID]. Based on the best-fit NFW profiles given in these works, the weak 
lensing measurements in [39] based on SDSS imaging and galaxy dynamics from [ID] 
give J values a factor of 2 higher while the weak lensing measurements from [38] based 
on CFHT imaging give J a factor of two lower than the value quoted in Table [U 

Based on the X-ray modeling, the nearby Fornax and M49 groups have the highest 
J values leading to the tightest constraints on dark matter models, so we specifically 
consider alternate determinations of the dark matter densities of these systems. The 
density profile of the central region of Fornax has been determined based on observations 
of the dynamics of globular clusters in its dominant central galaxy [UJ. Fig. |5] (left panel) 
compares the constraints on WIMP models for J as quoted in Table [T]for the smooth halo 
only, including the errors, to the corresponding constraints using J as determined from 
the best-fit NFW profile from the globular cluster sample in [41]. Here the dashed line 
shows the fit quoted in [41] for the full globular cluster sample while the dot-dashed line 
shows the fit to their "safe" subsample which employed a more conservative selection to 
reduce outliers; for these profiles, the expected gamma-ray flux from Fornax is reduced 
slightly by factors of 1.4 and 2.4, respectively. In Fig. (5j the comparison is only shown 
for a bb final state, but the same shift in normalization is expected for all annihilation 
final states. 

A determination of the dark matter density profile of M49 based on a combination 
of stellar kinematic and X-ray data is presented in [42] and compares quite well to the 
globular cluster kinematics presented in the same work. Fig. [5] (right panel) shows the 
comparison of the constraints for this profile. The density profile of [42J implies an 
enhancement of the gamma-ray flux relative to what is predicted from Table [I] by a 
factor of 2.3. 

In summary, determinations of J for some of the best candidate clusters in our 
sample based on different observables vary by roughly a factor of two, and the X-ray 
measurements we employ do not yield J values systematically higher or lower than 
other methods. In all cases, we have considered only the NFW functional form for the 
dark matter mass profile. However, since the intergral in Eq. ([2]) is taken over a large 
volume and for a very distant object, the assumed inner slope of the density profile 
does not have a large effect on the predicted gamma-ray flux. Alternate density profiles 
such as the cored Burkert profile [43] or the centrally steeper Moore profile [44] give 
predicted fluxes about 40% lower (Burkert) or higher (Moore) than the NFW profile, 
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Figure 5. Comparison of the annihilation cross-section constraints based on the X- 
ray modeling of the cluster density profiles described in Section 13.11 to the constraints 
using other observables for the two systems with the highest J, Fornax (left panel) and 
M49 (right panel). The solid lines and grey bands show the reference constraints and 
errors for J from Table [TJ For Fornax, these constraints are compared to modeling 
of the dark matter profile from globular cluster dynamics presented in [41) for both 
their full (dashed line) and "safe" (dot-dashed) globular cluster samples. For M49, 
we compare to the density profile determined by [42] through fitting a combination 
of stellar kinematic and X-ray data. In all cases, the dark matter density profile is 
parameterized with an NFW profile. A bb final state is assumed, and the dark matter 
models are the same as in the left panel of Fig. [3] 



similar to or smaller than the differences in flux presented above stemming from different 
observational determinations of the density profiles. 

4. Summary 

Clusters of galaxies are the most massive collapsed objects in the Universe and are very 
dark matter dominated, making them potentially bright sources of gamma-ray emission 
from dark matter annihilation. No significant gamma-ray emission has been detected 
from clusters of galaxies in the first 11 months of Fermi-LAT survey mode observations 
[10J. In this paper, we explored the implications of the non-detection of clusters by 
Fermi-LAT in terms of constraints on models of dark matter annihilation. In particular, 
we focused on the six best candidate clusters and groups of galaxies in the context of 
searches for gamma-ray emission from dark matter pair annihilation pj after excluding 
from the sample clusters which host bright central AGN or lie close to the Galactic plane. 
We analyzed the Fermi-LAT data to derive upper limits on the gamma-ray flux from 
dark matter annihilation in specific models, self-consistently incorporating the expected 
spectral shape for a given particle mass and annihilation final state. We conservatively 
assume only gamma-ray emission from dark matter annihilation when interpreting the 
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upper limits, though cosmic rays generated through astrophysical processes like shocks 
could also lead to gamma-rays from clusters. A much larger sample of clusters and a 
more general discussion of clusters as potential gamma-ray sources will be presented in 
a forthcoming paper. 

In particular, so-called lepto-philic dark matter models with relatively high particle 
masses and annihilation cross-sections proposed to fit the Pamela and Fermi-LAT 
e + e~ data would lead to significant IC gamma-ray emission in clusters [3]. Here we 
conservatively only consider IC scattering of the CMB. Including the contribution of 
substructure in clusters on galactic scales, known to exist, the non-detection of local 
clusters and groups excludes models fitting the Pamela positron excess with masses 
above 1 TeV for a jJ> + yr final state. Even for a smooth NFW dark matter profile with 
no substructure, models fitting the Pamela data with masses > 2 — 3 TeV are excluded. 
Unlike similar constraints placed by the non-detection of dwarf spheroidal galaxies by 
Fermi-LAT p3] which depend sensitively on the model employed for the spatial diffusion 
of e + e~, the leakage of cosmic rays due to spatial diffusion in clusters is not expected 
to be significant or significantly reduce the IC gamma-ray emission. 

With the inclusion of galactic scale substructure, the cluster gamma-ray constraints 
are beginning to probe some dark matter models annihilating to a bb or similar final state 
- producing a gamma-ray spectrum similar to what is expected in popular particle dark 
matter models, such as the neutralino - with thermal relic densities consistent with the 
observed universal matter density. Including the contribution from smaller mass dark 
matter substructures, the constraints we obtain become increasingly tighter. In general, 
however, these constraints are a bit weaker than (or in the most optimistic case similar 
to) the constraints derived from the non-detection by Fermi-LAT of the best candidate 
dwarf spheroidal galaxies [T5] . 

Clusters of galaxies, an as-yet gamma-ray quiet source class, will continue to 
be monitored as the Fermi Telescope increases its observing time in survey mode. 
In addition, future investigations will improve these constraints through stacking of 
clusters, improvements in our understanding of the gamma-ray backgrounds, and 
searches for extended emission around clusters like Perseus which host central gamma- 
ray AGN. While the detection of a gamma-ray signal from clusters might be attributed 
to a variety of source mechanisms, including cosmic rays from merger or accretion shocks 
and particle dark matter annihilation or decay, as shown in the present study even null 
results imply interesting and compelling constraints. 
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